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ABSTRACT 


An antenna engineering handbook and database program has been constructed by 
engineers at the Lawrence Livermore National Laboratory (LLNL) using the Numerical 
Electromagnetics Code (NEC) antenna modeling program to prepare data performance 
on tactical field communication antennas used by the Army. It is desirable to have this 
information installed on a personnel computer (PC), using relational database tech- 
niques to select antennas based on performance criteria. 

This thesis obtains and analyses current distributions and radiation pattern data by 
using NEC for the following set of four (4) high frequency (HF) tactical generic antennas 
to be used in future spherical mode expansion work: a quarter wavelength basic whip, 
a one-wavelength horizontal quad Loop, a 564-foot longwire, and a sloping “vee beam” 
dipole. 

The results of this study show that the basic whip antenna provides good ground 
waye communication, but it has poor near vertical incident skywave (NVIS) perform- 
ance. The current distribution has the characteristics of standing waves. The horizontal 
quad loop antenna is good for NVIS and medium range skywave communications. The 
current distribution is sinusoidal and continuous around the loop. The long wire an- 
tenna allows short, medium and long range communications and a standing wave cur- 
rent distribution occurs along the antenna axis due to non-termination. The sloping “vee 

beam” antenna favors long range communication and the current distribution is mainly 
that of travelling sinusoidal waves. --- ae < 

Because of their well-known efficiency, the basic whip and quad loop can be used 
as reference standards for the spherical mode expansion work. The longwire and sloping 


“vee beam” antenna are unwieldy, but they are effective as base station antennas. 
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THESIS DISCLAIMER 


This document was prepared as an account of work sponsored by an agency of the 
United States Government. Neither the United States Government nor anv of their 
employees, makes any warranty, express or implied, or assumes any legal hability or re- 
sponsibility for the accuracy, completeness, or usefulness of any information, apparatus, 
product, or process disclosed, or represents that its use would not infringe privately 
owned rights. Reference herein to anv specific commercial products, process, or service 
by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute 
or imply its endorsement, recommendation, or favoring by the United States Govern- 
ment. The views and opinions of authors expressed herein do not necessarily state or 
reflect those of the United States Government and shall not be used for advertising or 


product endorsement purposes. 
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Il. INTRODUCTION 


A. BACKGROUND 

A seven (7) volume antenna engineering handbook has been produced at the 
Lawrence Livermore National Laborotory (LLNL) using the Numerical 
Electromagnetics Code (NEC) antenna modeling program for the US Army Information 
Systems Engineering Center (USAISEC), Ft. Huachuca, Arizona. [Refs. 1,2]. 

A host of field tactical antenna types and configurations have been characterized for 
optimizing performance for high frequency (HF) communications when interfaced with 
existing propagation codes and models. 

It is desirable to have all useful antenna information installed in a database program 
on a personnel computer (PC) for selecting of best candidate antennas for a given op- 
erational scencrio. Initial work was started on this project several vears ago and has 
continued to date [Ref. 3]. 

Initial efforts here concentrated on applying existing software and new programs 
necessary for completing this effort and on producing a computer program for the PC. 

Most of the antennas in the antenna handbook are less than a few wavelengths in 
extent, limiting the directivity. This limitation translates into patterns which are 
describable in terms of a few low order spherical modes. 

A proposal, made by Dr. T. L. Simpson, Electrical and Computer Engineering De- 
partment, University of South Carolina, was to solve the problem of pattern generation 
by finding a more natural and simple way of expressing the currents in terms of a 
spherical modal expansion. “Since the general solution of the electromagnetic wave 
equation in spherical coordinates can be conveniently expressed as a sum over a Set of 
orthogonal eigenfunctions, or modes (Refs. 4, 5, 6], one can alwavs express a phvsically 
realizable pattern function such as an expansion. In effect, the use of the spatial filtering 
properties of the antenna reduces the antenna pattern to its fundamentally simplest 
form.” [Ref. 7]. 

It is desired to investigate and complete the representation of HF antenna patterns 
using Spherical Mode Expansion. The result of this work will be used to generate pat- 


terns on the PC with minimial disk storage requirements and fast execution speed. 


B. SCOPE OF THE THESIS 

In order to represent Hf{F antenna patterns using Spherical Mode Expansion, NEC 
results for radiation pattern calculations and current distribution data are required for 
the support of research in the spherical modal solutions. 

This thesis presents current distribution and radiation pattern data produced by 
NEC for the HF range of 5 and 30 MHz for the following set of four (4) tactical generic 
antennas: a quarter wavelength basic whip with 8’ ground stake, a one-wavelength hor- 


izontal quad loop, a 564-foot longwire, and a sloping “vee beam” dipole. 


C. SELECTION AND MODELING CRITERIA FOR THE ANTENNAS 
1. Ground and Frequency Selection 
Several options are available for modeling antennas which are located near the 
ground. Either perfect ground or imperfect ground can be specified, with the NEC user 
supplying values for the conductivity (a) and relative permittivity (¢,), in the case of finite 


ground, as shown by the following table: 


Table 1. GROUND TYPES 


parm 
Peon 5.0] stor and mein bile 










The ground constants (conductivity and relative permittivity) used throughout 
this thesis are “average” ground values characteristics of many rural sites, so that the 
performance calculated by the computer would be similar to the results one could expect 
to find under tvpical tactical conditions. 

The frequencies selected for this study are 5 MHz, which is the logarithmic 
mean for frequencies over which near vertical incident skywave (NVIS) communications 


are conducted, and 30 MHz, which is the highest frequency defined by the International 


Radio Consultative Commitee (CCIR) for the HF band (3 to 30 MHz). lonospheric 
skywave propagation is dominant at these frequencies and is used for skvwave links. 
The reason for selecting these frequencies is to calculate the performance of the 
generic antennas at the low and high frequency limits of the HF band. Thus the results 
bracket the performance over the entire HF band. 
2. Antenna Selection 
Several criteria were used when selecting the antennas for use in the Antenna 
Handbook. The following performance criteria were selected as desirable for field con- 
ditions: 
1. Communication effectiveness. 
2. Ruggedness. 


3. Easv of construction in the field. 


ae 


Simplicity. 


The Antenna Engineering Handbook presently contains the following set of 
antenna model configurations: 


1. Basic whip with and without a radio set. 
2. Basic half square. 

3. Inverted half square. 

4. Folded half square. 

5, 3-element half square. 

6. Inverted 3-element half square. 

7. Horizontal vee beam. 

§. Sloping vee beam. 

9. Sloping vee beam, legs terminated. 
10. Horizontal dipole. 

11. Inverted-vee dipole. 

12. Horizontal square loop (quad). 

13. Vertical triangular loop (delta Loop). 
14. Basic Jongwire. 

15. Sloping longwire. 


16. Inverted-vee longwire. 


From these, a set of four generic antennas mentioned in part B has been se- 
lected. 

Besides military use, some of the antennas selected for this work are commonly 
used by amateur, commercial, and maritime HF stations. These were selected because 
some of them have buried parts, others do not touch the ground, or the physical sizes 
differ substantially. Some of them are slanted while others have axes parallel to 
cartesian coordinates, etc. Thus the simulation program can be tested under different 
conditions for these well-known antennas. 


D. DESCRIPTION OF THE ANTENNAS 
1. Basic Whip 


e This antenna is a quarter wavelength long at the operating frequency. It is verti- 
callv polarized and is fed against ground. 


e The radius of the antenna is 0.015 m. 


e It has also an 8’ ground stake buried beneath the ground as seen in Figure I. 


2. One Wavelength Horizontal Quad Loop 


e The One wavelength quad loop antenna has one-quarter wavelength sides and is 
placed horizontally a quarter-wavelength high over the ground plane. The feed 
point is at one corner. 


¢ It has 0.001 m. radius. The structure is shown in Figure 2. 


3. 564-Foot Longwire 


e® This is a 564-foot long wire antenna with a 45° slanted feeder, and is placed 40 feet 
high above the ground. The feeder is driven by a voltage source whose lower end 
is tied to a 2-foot ground stake. 


e The radius of the wire is also 0.001 m. The structure is shown in Figure 3. 


4. Sloping “Vee Beam” Dipole 


e This antenna may be visualized as two diverging lines terminated through a 
600-ohm resistor to a ground stake. It is fed at the apex. 


e The apex height is 40 feet above the ground and both ends are 7 feet high. The 
length of each leg is 528 feet. 


¢ The conductor is also 0.001 m. in radius. The structure is shown in Figure 4. 


E. INTRODUCTION TO THE NUMERICAL ELECTROMAGNETICS CODE 
Computer modeling of antennas started in the late 1960's. Since then it has become 

a powerful tool for antenna design. Computer codes have been developed based on the 

Method-of-Moments, Geometrical Theory of Diffraction, or integration of Maxwell's 


QUARTER WAVE WHIP 


THETA = 60.00 PHI = 60.00 ETA = 90.00 





Figure 1. Quarter Wavelength Basic Whip with 8’ Ground Stake 





ONE WAVELENGTH HORIZONTAL QUAD LOOP 


QUARTER WAVELENGTH HIGH 


THETA = 60.00 PHi= 60.00 ETA = 90.00 





Figure 2. One Wavelength Horizontal Quad Loop 





564—FOOT LONGWIRE ANTENNA WITH SLANTED FEEDER 
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Figure 3. 564-Foot Longwire 

















SLOPING "'VEE BEAM" DIPOLE 


THETA = 60.00 PHI = 60.00 ETA = 90.00 


Figure 4. Sloping “Vee Beam” Dipole — 





Equations. The Numerical Electromagnetics code - Method of Moments [Ref. 2} has 
become the most widely used of all thin wire codes for modeling resonant-sized anten- 
nas. The reasons are that NEC features systematic updating and extension of its capa- 
bilities, user-oriented documentation and accessibility of its developers for user 
assistance. 

NEC has been under development for more than ten years (earlier known as 
BRACT and AMP). It is a hybrid code which uses an Electric Field Integral Equation 
(EFIE) to model wire like objects and a Magnetic Field Integral Equation (MFIE) to 
model closed surfaces with time harmonic excitation. NEC is commonly applied to 
modeling antennas in VLF-to-VHF applications on ships, vehicles or on the ground. 
A number of features are included for efficient modeling of antennas and scatterers in 
their environments, including antennas which are interacting with or buried in a finitely 
conducting ground. Recently, a capabilitv for modeling insulated wires in the air or 
ground has added to NEC and improved accuracy for electrically small antennas and the 
treatment of an abrupt change in wire radius [Ref. 8]. 

The utility of NEC has been enchanced by the development of an interactive 
graphics utility (IGUANA) [Ref. 9] for developing the model description and exanuning 
the results. IGUANA supports both NEC and the Method-of Moments Code 
MININEC [Ref. 10] for use on personnel computers. 








Il. PERFORMANCE RESULTS OF ANTENNAS 
A. GENERAL 


e All the antennas have been excited by a voltage source (applied E field source). 
e The antennas are driven by 1 kW input power. 
e The radiation pattern range is | km. 


e The wavelengths in the average ground have been calculated and found as 15.6 m 
for 5 MHz, and 3.1 m for 30 MHz. These data have been used to set segment 
lengths for buried ground stakes. 


e All of the calculations have been done in double precision at 32 bits for accuracy. 


B. BASIC WHIP WITH 8’ GROUND STAKE 

This antenna is basically a simple monopole with an 8’ ground stake. For spot fre- 
quency operation, a monopole antenna is satisfactory at a coastal site, but poor at low 
conductivity sites; its low angle coverage is limited by the amount of ground improve- 
ment. Therefore a ground stake has been added to the basic monopole structure. 

At 5 MHz, the antenna structure is 15 m above the ground ( 0.25 4 ), and is modeled 
by 6 segments, such that all the segment lengths are less than 6 m( 0.1 4 ). For the 
ground stake 3 segments were used due to the smaller wavelength in dirt. 

At 30 MHz, the length of the antenna itself is 2.5 m, and it has 5 segments. Also 
for the buried part. 15 segments wer used to model the ground stake. 

At 5 MHz, where the ground stake is electrically short, a part of the total base cur- 
rent comes from the return current which flows in the ground bevond the ground stake. 
The current distribution is shown in Figure 5. It is a standing wave current distribution. 
The current in the ground stake decrdases rapidly and flows outward in the ground. Also 
the phase is constant which is a characteristic of standing waves. 

For this current distribution, a radiation pattern has been calculated and plotted in 
Figure 6. The maximum field strength occurs at a 27° take-off angle. 


The radiation resistance was calculated as 





At 30 MHz the ground stake is about a quarter wavelength long. It provides good 
ground improvement favoring ground wave communication. The standing wave current 
distribution is shown in Figure 6. The phase difference between the whip antenna and 
the ground stake is 180° as expected. 

For this frequency, a radiation pattern has been calculated and is shown in Figure 
8. The maximum field strength occurs at 29° take-off angle. 

The radiation resistance at 30 MHz was found as 


ie 
Rie ae eee o (2) 


[7] (2.9 +j0.1)°| 


C. ONE WAVELENGTH HORIZONTAL QUAD LOOP 

The one wavelength quad loop antenna can be considered also as a modified folded 
dipole. By pulling the dipole wires apart at the center, the single-turn quad loop antenna 
is obtained. The length of each side is 0.254. 

At 5 MHz, the antenna is placed 15 m above the ground. Each side was modeled 
by 15 segments. At 30 MHz the antenna is 2.5 m above the ground and each side was 
modeled by 13 segments. 

For a one wavelength perimeter it is reasonable that the current distribution is 
sinusoidal as shown in Figure 9 for 5 MHz and in Figure 13 for 30 MHz. It is contin- 
uous around the loop and the magnitude is minimum at the corners between the feed 
point sides and the far sides. The phase reverses everv 0.5 /. 

Principal patterns at 5 MIIz are shown in the following order; Figure 10 shows the 
YZ-plane pattern plot of E,. In this plane the half-power beamwidth (IIPBW) is 76°; 
Figure 1] shows the XZ-plane pattern plot of E, and the HPBW is 118°; Figure 12 shows 
the NY-plane pattern plot of £,. 


The radiation resistunce was calculated as 


1000 


i TO LG: 3 
(esis e 


Pied = 

At 30 MHz the YZ-plane pattern plot of £, is shown in Figure 14 with HPBW = 

74°; Figure 15 shows the XZ-plane pattern plot of E, with HPBW = 116°; Figure 16 
shows the NXY-plane pattern plot of E,. The radiation resistance is 


1000 


OO = 132 
soso =oy (4) 
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QUARTER WAVE WHIP OVER AVE. GROUND / 5S MHZ 
WITH 8’ GROUND STAKE/INCLUDING THE GND STAKE CURRENT 


Current Magnitude 
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Figure 5. Current Distribution of Basic Whip Antenna at § MHz.: Bold line in- 


dicates the current distribution on the ground stake. 
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Figure 6. Elevation Pattern of Basic Whip Antenna at 5S MHz: 
€, = 10.0.6 = 0.003 








QUARTER WAVE WHIP WITH 8'GROUND STAKE / FREQ= 30 MHZ. 
OVER AVERAGE GROUND / INCLUDING THE GND. STAKE CURRENT 
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Figure 7. Current Distribution of Basic Whip Antenna at 30 MHz.: Bold line in- 


dicates the current distribution on the ground stake. 


QUARTER WAVE WHIP WITH 8' GROUND STAKE OVER AVERAGE GROUND. 


ELEVATION PATTERN / FREQ=30 MHZ. 
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ONE WAVELENGTH HORIZONTAL QUAD LOOP / FREQ= 5 MHZ. 
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Figure 9. Current Distribution of Quad Loop Antenna at 5 MHz. 
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Figure 10. The YZ-plane pattern of Quad Loop Antenna at 5 MUHz.: 
£,= 10.0, ¢ = 0.003 
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Figure 11. The XZ-plane pattern of Quad Loop Antenna at 5 MHz.: 
&,= 10.0, 0 = 0.003 








ONE WAVELENGTH HORIZONTAL QUAD LOOP / HEIGHT=.25 WAVELEN 
OVER AVERAGE GROUND / AZIMUTH PATTERN / FREQ=5 MHZ. 
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ONE WAVELENGTH HORIZONTAL QUAD LOOP / FREQ= 30 MHZ. 
OVER AVERAGE GROUND / HEIGHT= 0.25 WAVELENGTH. 


® 
2 
S 
im 
£2 
D 
a) 
= 
aie 
c 
® 
be 
pas 
3 
o 


25 50 75 100 125 180 175 200 225 250 275 
Electrical Distance (Degrees) 


- Current ‘Phase 


75 100 125 150 175 200 225 250 275 
Electrical Distance (Degrees) 





Figure 13. Current Distribution of Quad Loop Antenna at 30 MHz. 
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Figure 14. The YZ-plane pattern of Quad Loop Antenna at 30 MHz: 
&, = 10.0, 6 = 0.003 
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ONE WAVELENGTH HORIZONTAL QUAD LOOP / HEIGHT=.25 WAVELENGTH 
OVER AVERAGE GROUND / ELEVATION PATTERN PHI=0 / FREQ=30 MHZ. 





Figure 15. The XZ-plane pattern of Quad Loop Antenna at 30 MHz.: 
e,= 10.0, 6 = 0.003 
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D. LONG-WIRE ANTENNA 


The wire antennas introduced so far have been resonant-size antennas. The wave 


travelling outward and reflected from the end sets up a standing-wave type of current 
distribution. From this section on, non-resonant size antennas will be introduced. 

A long-wire antenna is one that is greater than one halfwave length long. When 
such antennas are unterminated at the end, the current distribution is mainly that of a 
standing wave. If the antenna is cut to a resonant length, the input impedance will be 
resistive. The theoretical current distribution is a sinusoidal standing wave. But because 
of the loss due to radiation, the actual radiation pattern differs from the theoretical cal- 
culation. The first recognizable difference is that the lobes tilt toward the unfed end. 
This difference between actual patterns and theoretical patterns (based on standing 
waves only) is much less in the case of center-fed antennas. 

If a longwire is terminated in a resistance equal to its characteristic impedance, the 
current distribution is mainly that of a travelling wave. The longer the wire, the smaller 
will be the angle between the wire and the main lobe. The important difference between 
the patterns of terminated and unterminated wire antennas is the absence of large rear 
lobes in the terminated case. 

In this study the long wire was 564 feet long and placed 40 feet above the ground. 
It was unterminated. 

At 5 MIiz the longwire is on the order of 3 wavelengths (171.9 m) and was modeled 
by 51 segments. Also this antenna has a 45° slanted feeder segmented by 8 pieces. This 
feeder was driven by a voltage source whose lower end is tied to a 2-foot ground stake. 
This ground stake was radius tapered and had 2 segments. At 30 MHz the antenna is 
approximately 17 wavelengths long and has 200 segments. The slanted feeder was 
modeled by 45 segments and the ground stake had 4 segments. 

For the unterminated case with a pure sinusoidal current distribution (standing 
wave), the source current can be expressed as a superposition of two uniform waves 


travelling in opposite directions as 


I(x) = I sin[mat ae + 0.5)] (5) 


where m Is an integer equal to Ak: L is the Jength of the antenna, and k is the 


Qn 
wavenumber (=). 
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Figure 17 shows the current distribution at 5 MHz. Standing waves occur along the 
antenna axis due to non-termination. Also the magnitude decreases slightly due to ra- 
diation. The phase also reverses every half wavelength. 

Figure 18 shows the elevation pattern at 5 MHz. It is symmetrical about the Z-axis 
and has a main beam calculated by 


Om = cos'(1 Sa) (6) 


which is 30°, and additional beams in every « = 180 — o,, direction providing good NVIS, 
medium range, and long range communications. 

Figure 19 shows the azimuth pattern for a 30° take-off angle. Figure 20 shows the 
same pattern for a 60° take-off angle. Figure 21 shows the azimuth pattern for an 80° 
take-off angle. All the radiation patterns above are symmetric along the antenna axes 
and the vertical gains are dominant to the horizontal gains. 

Radiation resistance at this frequency was calculated as 

S00 


Re 9474 (7) 
me 10.8 471.17 | 


Figure 22 shows the current distribution at 30 MHz. Standing waves exist as in the 
former case but the effect of losses is higher than for the same antenna at 5 MHz. In 
addition, there are also a few travelling waves. The minimum magnitude of the current 
occurs in the middle of the antenna. 

Figure 23 shows the elevation pattern at 30 MHz. It has a main beam at 13°. Al- 
though the radiation pattern has symmetry about the z-axis, the gain of the front lobe 
is larger than the back lobe. 

Figure 24 and 25 show the azimuth patterns for 15° and 25° take-off angles respec- 
uvely. Both radiation patterns show the main lobe in the forward direction larger than 
in the reverse direction. This indicates also the existence of travelling waves on the an- 
tenna. 


Radiation resistance at this frequency was calculated as 


500 


ILI &. (8) 
(1.1 — 1.2)" 


rad = 








564—FOOT LONGWIRE WITH SLANTED FEEDER / FREQ= 5 MHZ. 
OVER AVERAGE GROUND / HEIGHT = 40 FEET. 
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Figure 17. Current Distribution of Long Wire Antenna at 5 MHz. 
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564—FOOT LONGWIRE WITH SLANTED FEEDER / HEIGHT=40 FEET. 
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Figure 18. Elevation Pattern of Long Wire Antenna at §S MHz: 
&, = 10.0, ¢ = 0.003 
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564—FOOT LONGWIRE WITH SLANTED FEEDER / FREQ= 30 MHZ. 
OVER AVERAGE GROUND / HEIGHT = 40 FEET. 
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Figure 22. Current Distribution of Long Wire Antenna at 30 MHz. 
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Figure 23. Elevation Pattern of Long Wire Antenna at 30 MHz.: 
&, = 10.0, 6 = 0.003 
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E. SLOPING VEE BEAM ANTENNA 

This is a travelling-wave, balanced antenna. The angle between the descending wires 
was chosen so that the principal lobes of two radiators are superimposed. The mast 
height is chosen to make the antenna tilt angle the same as the take-off angle. Half of 
the radiated power is reflected from the ground. Therefore the ground should be smooth 
and free from obstructions. This antenna is horizontally polarized and has a relatively 
high input impedance. 

The modeled antenna had two 528-foot descending wires. The interior angle be- 
tween two wires was 60°. The feed point was 40 feet above the ground and each leg of 
the “Vee” was terminated through a 600 Q resistor to a 2 feet ground stake. 

At 5 MHz each 528-foot leg has 45 segments. The load is placed on one segment 
at this frequency. The ground stake was modeled by 2 segments. 

At 30 MHz each leg had 221 segments. The wire connecting the load to the ground 
stake has 4 tapered-length segments. Each ground stake was modeled with 4 segments. 
The gain at lower take-off angles provides good long range communication. Figure 33 
shows the azimuth pattern at a 10° take-off angle. Although it is horizontally polarized, 
some of the radiated power comes from ground reflection. 

For the current distribution on antenna wire No. 1, 


1, =1,e°" (9) 
and for the current distribution on antenna wire No. 2, 
Lh=1,e" 


Where /, is the current amplitude and s is the distance from the feed point. 

Figure 26 shows that the current distribution at 5 MHz is essentially a travelling 
wave sinusoidal current distribution. The magnitude decreases gradually toward the 
load and the phase reverses every half wavelength. 

Figure 27 shows the elevation pattern of E, at @ = 60°, the lobe created by one 
descending wire. Figure 28 shows the YZ-plane pattern plot of E,. The lobe nearest to 
the antenna wire is the largest (the main lobe) and all others progressively diminish in 
strength. Figure 29 shows the azimuth plane pattern, £,, at an elevation angle of 20° . 
The horizontal pattern gain in the direction of antenna is maximum and the back lobes 
are introduced by ground reflections. The input impedance is highly resistive. At this 


frequency the radiation resistance was calculated as 


Los) 
‘an 


1000 


(1 + 70.006) —1.1—Jo.006), ~ 8? @ (10) 


Riad = 

Figure 30 demonstrates that the current distribution at 30 MHz is also a sinusoidal 
distribution with logarithmically decreasing magnitude. 

At 30 MHz the antenna is about 16 wavelengths long. Figure 31 shows the ele- 
vation pattern at @ = 60°. Figure 32 shows also the YZ-plane pattern of E,. Here the 
principles lobes are superimposed. The gain is maximum at a low take-off angle pro- 
viding good long range communication. Figure 33 shows the azimuth pattern at a 10° 
take-off angle. 


At this frequency radiation resistance was calculated as 


1000 
Rraa = (1.2 +j0.005)( —1.2 —j0.005)|_~ ae ay) 
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Figure 26. Current Distribution of Sloping Vee Antenna at 5 MHz. 
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SLOPING "VEE BEAM" DIPOLE / LENGTH OF EACH LEG = 528 FEET 
OVER AVERAGE GROUND / ELEVATION PATTERN PHI=60 / FREQ=5 MHZ. 
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Figure 27. Elevation Pattern of Sloping Vee Antenna at 5 MHz. ¢ = 60°: 
e, = 10.0, ¢ = 0.003 
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Figure 28. YZ-Plane Pattern of Sloping Vee Antenna at 5 MHz: 
é, = 10.0, 6 = 0.003 
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SLOPING ''VEE BEAM" DIPOLE / ONE LEG / FREQ=30 MHZ. 
OVER AVERAGE GROUND / INCLUDING THE GND. STAKE CURRENT 


Current Magnitude 


0 1090 2000 3000 4000 5000 6000 
Electrical Distance (Degrees) 


200 


100 


Current Phase 
o 


) 1000 2000 3000 4000 5000 6000 
Electrical Distance (Degrees) 


Figure 30. Current Distribution of Sloping Vee Antenna at 30 MHz. 
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Figure 31. Elevation Pattern of Sloping Vee Antenna at 30 MHz. ¢ = 60°: 
e, = 10.0, 6 = 0.003 
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_ Figure 32. YZ-Plane Pattern of Sloping Vee Antenna at 30 MHz: 
e, = 10.0,¢ = 0.003 
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I. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This thesis investigated the performance of 4 tactical HF field antennas at 5 and 30 
MHz over average ground conditions. 

The basic whip antenna provides good ground wave communication. The gain of 
the ground wave is proportional to the ground conductivity improvement. The skywave 
is not as dependant on ground constants. It has poor NVIS performance, but its low 
take-off angle favors long distance communication. At HF the usefulness of the basic 
whip is limited by the wide operating bandwidth usually required. The current distrib- 
ution has the characteristics of standing waves. 

The horizontal quad loop antenna has highest gain at high take-off angles. There- 
fore it has very poor ground wave performance. The radiation patterns show that it is 
good for NVIS and medium range skywave communications. The one wavelength quad 
loop has the same elevation and azimuth pattern as a half wavelength dipole, but it has 
more gain. Also the loop antenna has lower terminal resistance than a folded dipole. 
The current distribution is sinusoidal and continuous around the loop. 

The long wire antenna allows short, medium, and long range communications. The 
radiation patterns at different frequencies show that as the frequency increases, the main 
lobe becomes more lobey and has better vertical gain at lower angles, providing long- 
range communication. The performance of the antenna remains practically unchanged 
within 10 MHz of the design center frequency. Generally, antennas which use termi- 
nated lines appear to have more radiation resistance than unterminated ones. Radiation 
and wire losses and the presence of ground and obstacles affect the pattern of anv ideal 
source. If improvement is required, it is better to lengthen the longwire than to raise it. 
The height of the antenna can be adjusted so that the reflected wire combines in phase 
with the direct wave in the direction of maximum radiation. The important consider- 
ation here is the phase shift at the point of reflection [Ref. 11]. Since it is unterminated, 
a standing wave current distribution occurs along the antenna axis. 

For the sloping “vee beam” antenna, it is also true that, like a longwire, this antenna 
has also a relatively wide bandwidth. Since half of the radiated power reflects from 
ground, the ground should be free fr.n obstacles. The level of the sidelobes change 


considerably, depending on the type of ground. The level of the largest sidelobes is only 


f. 
SN 


about 4 - 7 dB below the main beam. This is why the usefulness of the antenna is re- 
stricted. Its low take-off angle provides very good long-range communication. The 
current distribution is mainly that of travelling sinusoidal waves. 

As a result, because of their well-known efficiency, the basic whip and quad loop 
can be used as reference standards for the spherical mode expansion work. The longwire 
and 528-foot sloping “vee beam ” dipole are large, unwieldy and difficult to construct, 
but they are effective as base station antennas. 


B. RECOMMENDATIONS 

The configurations of different antennas can be studied under different conditions 
or for different frequencies to provide more information for spherical modal expansion 
work, if needed. 
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APPENDIX A. NEC AND AUXILIARY PROGRAMS 


A brief description of some of the programs found in the NEC FORTRAN library 
are described below: 
1. SOMNTX FOR SOMMERFELD / NORTON GROUND METHOD 
This option can be requested on the GN card. This method is available for wires 
only. It uses an exact solution for the fields in the presence of ground and 1s accurate 
close to and beneath the ground. NEC reads interpolation tables from the file TAPE21. 
This input file must be created by running SOMNTX. SOMNTX reads 4 parameters 
from a data file: 
¢ EPR = relative dielectric constant of ground ( ¢, ) 
© SIG = conductivity of ground in mho;m(¢ ) 
© FMHZ = frequency in MHz. 


e¢ IPT = print option of interpolation tables. 


2. THE NUMERICAL GREEN’S FUNCTION (NGF) OPTION 
The main purpose of the NGF is to avoid unnecessary repetition of calculations. 
Another reason for using the NGF option is to exploit partial symmetry in a structure. 
This option can be requested by a WG card and a GF card. 
3. DNPGNEC 
A double precison version of NEC-? available in mainframe at NPS, it is good 
for modeling structures up to 300 segmem.s. DNPG1000, another version of NEC, is 
capable of handling up to 1000 segments. 
4. PLOT PROGRAMS 
¢ PLOTDGLP -- Plots data geometry to a laser printer. 
¢ PLOTDG79 -- Plots data geometry to a 3279 terminal. 


¢ PLOTRPHM -- Plots horizontal and vertical components of radiation patterrs. 
Use this to get azimuth plots in the horizontal plane. 


¢ PLOTRPVM -- Plots radiation pattern in the vertical plane. 


¢ CURRPLOT -- This program plots current magnitute and phase for selected seg- 
ments and tags or selected groups of segments or tags. 
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5. RVAL 
It is used for variable segmentation or tapered radius or both. RVAL asks the 
length of the first and the last segments and the total length of the wire. Then it calcu- 
lates the segment lengths and the segmentation ratio. It is used on GC cards, which are 
read only when the radius on a GW card is set to zero. 
6. DIRTLAM 
This program calculates the wavelength in dirt. The result is used to set segment 
lengths for wires buried under ground. It asks the relative dielectric constant, 
conductivity and operation frequency and gives the wavelength in the ground. 
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APPENDIX B. NEC DATA SETS FOR CURRENTS 


GEOMETRY : QUARTER WAVELENGTH BASIC WHIP WITH 8' GROUND ROD 
RADIO : NONE 

FREQUENCY : 5 MEGAHERTZ 

GROUND : EPSILON = 10; SIGMA = 0.003 


»9,0,0,0,0,14. 989,0. 015 
»0,0,0,0,0,-2. 438,0. 015 


GEOMETRY : QUARTER WAVELENGTH BASIC WHIP WITH 8' GROUND ROD 
RADIO : NONE 

FREQUENCY : 30 MEGAHERTZ 

GROUND : EPSILON = 10; SIGMA = .003 


0,0,0, 0,0,2.498, 0.015 
50,0, 0,0,-2.438, 0.015 


GEOMETRY : ONE WAVELENGTH HORIZONTAL QUAD LOOP 
HEIGHT = 0.25 WAVELENGTH 
CONDUCTOR = #12 AWG 

FREQUENCY : 5 MEGAHERTZ 

GROUND : EPSILON = 10; SIGMA = .003 


1,12, 0,0,14.989, -10.599,10.599,14.989, 0.001 
2,12, -10.599,10.599,14.989, 0,21.198,14.989, 0.001 
3,12, 0,0,14.989, 10.599,10.599,14.989, 0.001 

4,12, 10.599,10.599,14.989, 0,21.198,14.989, 0.001 


0,0,0,0, 5, 
2,0,0,0, 10,.003 
1,2,0,0 

0,1,1,01, 251.115 
0,3,1,01, -251.115 
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GEOMETRY : ONE WAVELENGTH HORIZONTAL QUAD LOOP 
HEIGHT = 0.25 WAVELENGTH 
CONDUCTOR = #12 AWG 

FREQUENCY : 30 MEGAHERTZ 

GROUND : EPSILON = 10; SIGMA = .003 


1,13, 0,0,2.498, -1.766,1. 766,2.498, 0.001 
2,13, -1. 766,1. 766,2.498, 0,3. 532,2.498, 0.001 
3,13, 0,0,2.498, 1. 766,1. 766,2.498, 0.001 

4,13, 1.766,1. 766,2.498, 0,3.532,2.498, 0.001 


2,0,0,0, 10,.003 

i 2 0,0 

0,1,1,01, 252. 345 

0,3,1,01, -252. 345 

GEOMETRY : 564-FOOT LONGWIRE WITH 45 DEGREE SLANTED FEEDER 
HEIGHT ABOVE GROUND = 40 FEET 
CONDUCTOR IS # 12 AWG 
FEEDER IS DRIVEN BY A VOLTAGE SOURCE WHOSE LOWER 
END IS TIED TO A 2-FOOT GROUND STAKE 

FREQUENCY : 5 MHZ 

GROUND =: EPSILON = 10 ; SIGMA = .003 

1,2, 0.0762,0. 06096,-0. 6096, 0.0762,0. 06096,0 

0,0, 1, 0.004,0. 002 

2,1, 0.0762,0. 06096,0, 0.0762,0.06096,1.302, 0.001 

3,8, 0.0762,0. 06096,1. 302, -10. 8138,0. 06096,12.192, 0.001 

4,51, -10. 8138,0. 06096,12.192, -182.721,0.06096,12.192, 0.001 

-1,2 

0,0,0,0, 5 

2,0,0,0, 10,.002 

1,2,0,0 

0,3,1,01, 1204. 350 

GEOMETRY : 564-FOOT LONGWIRE WITH 45 DEGREE SLANTED FEEDER 
HEIGHT ABOVE GROUND = 40 FEET 
CONDUCTOR IS # 12 AWG 
FEEDER IS DRIVEN BY A VOLTAGE SOURCE WHOSE LOWER 
END IS TIED TO A 2-FOOT GROUND STAKE 

FREQUENCY : 30 MHZ 


GROUND : EPSILON = 10 ; SIGMA = .003 


1,4, 0.0762,0. 06096,-0. 6096, 0.0762,0. 06096,0 











0,0, 1, 0.004,0. 002 
2,6, 0.0762,0.06096,0, 0.0762,0. 06096,1. 302, 0.001 
3,45, 0.0762,0.06096,1. 302, -10. 8138,0. 06096,12.192, 0.001 
4,200, -10. 8138,0. 06096,12.192, -182.721,0.06096,12.192, 0.001 
-1,2 
0,0,0,0, 30 
2,0,0,0, 19,.003 
1,2,0,0 
0,3,1,01, 872.360 
GEOMETRY : SLOPING "VEE BEAM" DIPOLE 
APEX HEIGHT = 40 FEET 
BOTH ENDS AT HEIGHT = 7 FEET 
LENGTH OF EACH LEG OF "VEE" = 528 FEET 
EACH LEG OF "VEE" IS TERMINATED THROUGH A 600-OHM 
RESISTOR TO A 2-FOOT GROUND STAKE 
INTERIOR ANGLE OF "VEE" = 60 DEGREES 
CONDUCTOR IS # 12 AWG 
FREQUENCY: 5 MHZ 
GROUND : EPSILON = 10 ; SIGMA = .003 
1,45, 0,0,12.192, 80. 3099,139. 1008,2. 1336, 0.001 
2,1, 80. 3099,139. 1008,2. 1336, 80. 3099,139.1008,0, 0.001 
3,2, 80. 3099,139.1008,0, 80. 3099,139. 1008, -0. 6096 
0,0, 1, 0.002,0. 004 
4,45, 0,0,12.192, -80. 3099,139. 1008,2. 1336, 0.001 
5,1, -80. 3099,139. 1008,2. 1336, -80.3099,139.1008,0, 0.001 
6,2, -80. 3099,139.1008,0, -80. 3099,139. 1008, -0. 6096 
0,0, 1, 0.002,0. 004 
122 
4, 2,1,1, 600 
4, 5,1,1, 600 
0,0,0,0, 5 
2,0,0,0, 10,.003 
1,2,0,0 
0,1,1,01, 436.485 
0,4,1,01, -436. 485 
GEOMETRY : SLOPING "VEE BEAM" DIPOLE 
APEX HEIGHT = 40 FEET 
BOTH ENDS AT HEIGHT = 7 FEET 
LENGTH OF EACH LEG OF "VEE" = 528 FEET 
EACH LEG OF "VEE" IS TERMINATED THROUGH A 600-OHM 
RESISTOR TO A 2-FOOT GROUND STAKE 
INTERIOR ANGLE OF "VEE" = 60 DEGREES 
CONDUCTOR IS # 12 AWG 
FREQUENCY: 5 MHZ 
GROUND : EPSILON = 10 ; SIGMA = .003 
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1,221, 0,0,12.192, 80. 3099,139. 1008,2. 1336, 0.001 
2,4, 80. 3099,139. 1008,2. 1336, 80. 3099,139. 1008,0, 
0,0,1. 306,0. 001,0. 001 

3,4, 80. 3099,139. 1008,0, 80. 3099,139. 1008,-0. 6096 

, 1, 0.002,0. 004 

21, 0,0,12.192, -80. 3099,139. 1008,2. 1336, 0.001 
-80. 3099,139. 1008,2. 1336, -80. 3099,139. 1008,0, 
1. 306,0. 001,0. 001 

~80. 3099,139. 1008,0, -80. 3099,139. 1008, -0. 6096 
1, 0.002,0. 004 


1, 415.428 
1 


APPENDIX C. NUMERICAL METHODS AND TECHNIQUES FOR 
ANTENNA ANALYSIS 


A. THE METHOD OF SOLUTION 

NEC uses both an electric-field integral equation (EFIE) and a magnetic-field inte- 
gral equation (MFIE). The EFIE in NEC is well-suited for thin-wire structures of small 
or vanishing conductor volume. MFIE is only useful for closed structures, especially 
those having large smooth surfaces. For a structure containing both wires and surfaces, 
the EFIE and MFIE are coupled in NEC. 

To model wires, the thin Wire approximation is used in the EFIE to reduce it to a 
scalar integral equation. The assumptions used are [Ref 8]: 


1. Transverse currents are neglected relative to axial currents on the wire. 


tr 


The circumferential variation in the axial current is neglected. 


Lo>) 


. The current is represented by a filament on the wire axis. 


4. The boundary condition on the electric field is enforced in the axial direction only. 


These approximations are valid for the wire radius much less than the wavelength 
and much less than the wire Iength. An alternative kernel for the EFIE, based on an 
extended thin-wire approximation, is also included in NEC for wires having too large a 
radius for the thin-wire approximation [Ref. 12]. With the thin-wire approximation the 
EFIE becomes: 


AN AN 


=. E()=——] 168 Re == 
: cr ce 








<— )g(F, 7") ds’ (12) 


where 


= unit vector along the wire axis. 


distance along the wire axis. 
E'(F) = incident electric field at r. 
@ = 2zf, where fis frequency. 


axial current. 


_ 

“, 

— 
It 


ll 


permittivity. 
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k = w,/ue = phase constant. 


observation point. 
r’ source point. 

g(7,r’) = exp = = free space Green’s function. 
R = (7—-7r’') 


The MFIE for a closed surface § is: 


Jf) = 2n x HF) + sth x [20 x V'gds'res' 


5 


JM?) surface current density. 


H(?) = incident magnetic field at the observation point. 


n unit vector normal to the surface at 7. 

The integral equations (12) and (13) are solved numerically in NEC by a form of the 
method of moments. This solution includes expanding the unknown current in a sum- 
mation of basis functions with unknown amplitudes and enforcing equality of weighted 
integrals of fields. The weighting functions for both wires and surfaces in NEC are 
chosen to be delta functions which resulted in a point sampling of the fields. Wires are 
divided into short segments with a sample point at the center of each segment. The sur- 
faces are approximated by a set of flat patches or facets with a sample point at the center 
of each patch. 

Delta functions are also used in NEC as the current expansion functions on a sur- 
faces. For the MFIE, Galerkin’s method [Ref. 13] provides good accuracy on large 
smooth surfaces. Due to nature of the integral-equation kernels, however, the choice 
of current expansion functions is more critical in the EFIE than the MFIE. 

A more realistic representation of the surface current is needed where a wire con- 
nects to a surface than the delta function expansion normally used with the MFIE. In 
the region of the wire connection, the surface current includes a singular component due 
to the current flowing from the wire onto the surface. 

As a.result of the continuity conditions on current and charge there remains one 
unknown associated with each basis function to be determined by solving the 


electromagnetic interaction equations. Enforcing appropriate conditions on the current 
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and charge results in rapid convergence of the solution. The conditions may be inap- 
propriate in some cases, such as for short segments at a change in radius and at junc- 
tions of closely spaced wires. The result is then slower convergence. More accurate 
conditions for these cases are currently being developed for NEC [Ref. 8]. 


B. CAPABILITIES OF NEC 
There are basically two versions of NEC written in FORTRAN. NEC-2 has no 
distribution restrictions and limited a maximum of 300 segments to model a given 
structure. On the other hand NEC-3 is classified as Defense Critical Technology and 
can be released only to DoD agencies and their contractors without specific DoD ap- 
proval. NEC-3 can be used up to 300, 1000, or 2000 segments, depending on the ver- 
sions of NEC being used. The capabilities of these codes are summarized below. 
1. Input 
A user-oriented data card set for a single run consists of three types of data 
cards. The deck must begin with one or more comment cards which can contain a brief 
description anc structure parameters for the run. These are followed by geometry data 
cards which permits defining straight wires, arcs, and surfaces. Electrical connections 
are determined in the program by searching for wire ends and patch centers that coin- 
cide. Shifts, rotation, and reflections can be used in building structures. Finally, a sec- 
tion of program control cards specifies electrical parameters such as frequency, loading 
and excitation, and requests calculation of antenna currents and fields. 
2. Output 
Outputs are selectable by input parameters and may include one or more fol- 
lowing data: 
1. charge density on wires. 
2. current. 
3. input impedance, admittance, and power. 
4. radiated power, ohmic loss, efficiency. 
5. radiated field, power gain and directive gain. 
average gain. 
near E and H fields. 


maximum coupling. 


eo Rr 


receiving patterns. 


10. scattering cross sections. 
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3. Source modeling 
A structure may be excited by one of the following sources: 


1. voltage source (applied-E-field source). 

. incident plane wave, linear polarization. 

. incident plane wave, right-hand elliptic polarization. 
. incident plane wave, left-hand elliptic polarization. 


. elementary current source. 


Roe TA BB WwW WY 


. Voltage source (current-slope-discontinuity). 


4. Nonradiating Networks and Transmission lines 
Nonradiating two-port networks and transmission lines may connect points on 
wires. The defining parameters are characteristic impedance, length, and shunt 
admittance. 
5. Loading 
Lumped or distributed RLC loads may be specified on wires. The options that 
the user can specif¥ are: 
1. series RLC, input ohms, henries, farads. 
2. parellel RLC, input ohms, henries, farads. 


series RLC, input ohms meter, henries meter, farads/meter. 


bow 


parallel RLC, input ohms meter, henries meter, farads meter. 


Also the conductivity of a round wire may be specified and the impedance 

computed taking account of skin depth. 
6. Ground Effects 

For an antenna structure in or near the ground, there are three options avail- 
able. A perfectly conducting ground is modeled by including the image field in the kernel 
of the integral equations. An approximate model for a finitely conducting ground uses 
the image modified by Fresnel reflection coefficients. This approximation is usable for 
antennas at least 0.1 to 0.2 wavelengths above the ground. 

NEC-3 includes an accurate treatment for wire structures above, below, or 
penetrating the ground surface. The solution is based on the Sommerfeld integral for- 


mulation for the field near the interface. 


7. Modeling Guidelines 

Some of the guidelines for modeling are summarized below. 

For accurate results, the lengths of wire segments should be less than about 0.1 
4. Longer segments up to about 0.14 4 may be acceptable on long straight wires or 
noncritical parts of a structure. For critical regions, shorter segments 0.05 / or less may 
be needed. The minimum segment length set by the computational precision is about 
10-44. In double precision the minimum length is about 10-*/. 

The wire radius, a, relative to / is limited by approximations used in the kernel 
of the EFIE. Two options are available, the thin-wire kernel and the extended thin-wire 
kerne!. With the extended thin-wire Kernel the ratio of segment length to radius can be 
as small as 0.5. With the thin-wire approximation, the limit is 2. In either case, only 
currents in the axial direction on a segment are considered, and there is no allowance for 
Variation of the current around the wire circumference. The acceptability of these ap- 
proximations depends on both the value of a4 and the tendency of the excitation to 
produce circumferential currents. Discontinuites in wire radius should be kept to less 
than a factor of two. 

Electricaily small loops present special accuracy problems. In single precision 
the minimum loop circumference is about 0.002 4 ; while in double precision it is about 
10-4, , 

The surfaces modeled with patches must be closed. A reasonable maximum for 
the area of a surface patch appears to be 0.04 square wavelengths. 

An important consideration in using NEC is the solution time versus model size 
since this may limit the amount of detail that can be modeled and the segment and patch 
densities. For a model using N wire segments, the solution time in seconds on a CDC 


computer is approximately [Ref. 8 ] 
T = 310 )AN7/ AM + UO NSLP (14) 


where M is the number of degrees of symmetry and k is | for free space, 2 for perfectly 
conducting ground modeled with the reflection coefficient approximation, and 4 to § for 
Sommerfield integral treatment of finitely conducting ground. 

The first term in Eq. (3) is the time to fill the interaction matrix and the second 
term is the time to factor it into triangular parts for solution. For a model using \, 


patches, the solution time is about [Ref. 8] 


T = (10™°)A(2.N,Y/M + 200° )(2.N,) M7 (15) 





since two rows and columns in the matrix are associated with each patch. 

When the Sommerfeld ground treatment is used, a fixed time of about 15 sec- 
onds [Ref. 8] is needed to generate the interpolation tables. The tables depend only on 
the ground parameters and frequency, and can be saved for use in any case in which 


these parameters are the same. 


C. CONCLUSION 

NEC is a versatile code for modeling antennas and their environment including 
transmission lines, networks, loading and ground effects. It is written in FORTRAN 
and consists of about 10,000 lines of code. 

Recently released extensions of NEC include a model for insulated wires and im- 
proved accuracy in modeling electrically small antennas and discontinuities in wire ra- 
dius. A long-range goal is to incorporate expert-systems technology to simplify the 
process of building models that conform to the modeling guidelines and to verify that 
the code is producing correct results {Ref. 8]. 

NEC will continue to be a widely used antenna modeling program, due to its con- 


venient operation, documentation, availability, and continuing support. 
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